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Pulsed  Thruster  Thrust  Stand  Measurement  Evaluations 

Keith  A.  McFalf ,  Gregory  G.  Spanjers\  John  H.  Schilling* 

Propulsion  Directorate,  Air  Force  Research  Laboratory 
Edwards  AFB,  CA  93524 

This  paper  presents  a  numerical  analysis  of  thrust  stand  response  during  pulsed  thruster 
operation.  The  system  model  incorporates  an  underdamped  displacement-type  thrust  stand  under  a 
periodic,  impulsive  load.  Analytic  and  numerical  methods  are  used  to  examine  thrust  stan 
measurement  characteristics  during  simulated  operation.  Displacement  oscillation  characteristics  as 
a  function  of  damping  coefficient  (a)  and  the  ratio  (t)  between  the  time  between  thruster  impulses 
and  thrust  stand  natural  period  are  presented.  The  effects  of  simulated  thrust  stand  random 
displacement  signal  noise  were  quantified  for  <x=0.3  with  X  varying  from  0.1  to  0.7  These  resu  ts 
show  that  measurement  uncertainty  due  to  noise  decreases  when  xis  increased.  An  analytic 
expression  relating  average  thrust  and  average  thrust  stand  displacement  was  derived.  Using  this 
expression,  the  effect  of  time-dependent,  random  thruster  impulse  variability  was  examine  . 
Numerical  simulations  were  used  to  estimate  the  dependence  of  thrust  measurement  uncertainty  on 
the  number  of  impulses  used  for  averaging.  These  results  also  showed  that  measurement  uncertainty 
due  to  impulse  bit  variability  decreased  when  x  is  increased.  Integration  of  10  impulses  is  adequate  to 
achieve  measurement  uncertainty  below  2%  for  a  random  10%  impulse  bit  variation. 
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NOMENCLATURE 

=  displacement  coefficients 
=  displacement  coefficients  for  uniform 
impulses 

=  fractional  displacement  coefficient 
differences:  equal  to  [ck  -Ck  Jl  Ck 
=  fractional  time  within  an  impulse 
=  thruster  impulse  bit 
=  index  for  counting  impulses 
=  thrust  stand  spring  constant 
=  thrust  stand  effective  mass 
=  number  of  impulses  averaged 
=  time 

-  thruster  pulse  period:  time  between  impulses 
=  natural  period:  equal  to  2 K  /  »JK  hn 

-  displacement  amplitude 
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time  averaged  displacement 
maximum  displacement  during  impulse 
displacement  at  time  of  impulse 
dimensionless  linear  independent  solution  to 
the  thrust  stand  equation  of  motion 
damping  coefficient  and  damping  rate 
fractional  displacement  noise 
fractional  impulse  bit  variability 
noise  and  impulse  induced  measurement 
uncertainty 

amplitude  fraction  corresponding  to  average 
thrust 

pulse  period  ratio:  equal  to  tp  /  Tnat 


INTRODUCTION 

The  Air  Force  Research  Laboratory  is 
investigating  the  application  of  small  satellite  propulsion 
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technology  to  support  future  DOD  propulsion 
requirements1'2.  NASA  has  evaluated  the  use  of  pulsed 
plasma  thrusters  (PPT)  for  the  proposed  third  deep  space 
mission  (DS-3)  of  the  New  Millenium  Program3.  With  the 
continuing  interest  in  these  and  other  small  satellite 
missions4  5,  propulsion  system  designs  are  being  explored 
which  use  pulsed  operation  to  achieve  both  high  specific 
impulse  and  low  power  operation.  An  accurate 
measurement  of  thruster  performance  is  critical  for 

mission  analysis  and  thruster  selection. 

A  pulsed  thruster  applying  an  impulsive  load  on 
a  displacement-type  thrust  stand6'7  8  results  in  oscillations 
about  the  average  displacement.  When  the  thruster 
discharge  frequency  is  much  greater  than  the  natural 
frequency  of  the  thrust  stand,  and  the  discharges  produce 
uniform  (shot  to  shot)  impulse  bits,  the  thrust  stand 
displacement  is  nearly  fixed  as  a  function  of  time. 
However,  as  the  discharge  frequency  is  reduced,  the 
amplitude  of  the  displacement  oscillation  increases. 

Methods  to  accurately  measure  the  impulse  bit 
have  included  rapidly  pulsed  thruster  operation  (where  the 
ratio  (t)  between  the  time  between  successive  impulses 
(thruster  pulse  period,  tr)  and  thrust  stand  natural  period 
(Tnm)  is  very  small,  *  1  /  30)6,  the  use  of  single  pulse 
impulse  bit  measurements6'78,  or  the  use  of  very  long 
period  thrust  stand  designs8.  These  devices  have  measured 
impulse  bits  in  the  100-300  pN-s  range,  with 
measurement  uncertainties  ranging  from  ±5%  to  less 
than  2%7.  While  these  devices  have  demonstrated  a  high 
degree  of  measurement  accuracy,  the  analytical  models 
used  to  interpret  the  data  are  only  strictly  valid  for 
measurements  with  r  «  1  or  with  a  single  impulse.  Past 
modeling  efforts  have  examined  the  transformation  of  the 
equations  of  motion  for  a  torsional  thrust  stand  from  a 
rotational  to  a  translational  reference  frame6,  damped  and 
undamped9  single  impulse  displacement,  and  the  very 
long  period  motion  of  “Watt  s  Pendulum. 

This  paper  presents  an  analysis  of  the  dynamic 
response  of  an  under-damped  displacement-type  thrust 
stand  under  periodic,  impulsive  loads.  Analytic  and 
numerical  methods  are  used  to  examine  thrust  stand 
displacement  characteristics  during  simulated  operation 
for  a  range  of  damping  coefficients  (a)  and  t.  Thrust 
measurement  approximations  and  accuracy  are  presented. 
The  effects  of  simulated  thrust  stand  displacement  noise 
and  time-dependent  random  thruster  impulse  variability 
are  examined. 

The  application  of  a  comprehensive  theory  of 
thrust  stand  operation  under  periodic,  impulsive  loads  has 
the  potential  to  provide  accurate  measurements  over  a 
wide  range  of  X.  Such  a  capability  would  allow  greater 
flexibility  in  experiment  design.  In  addition,  an  improved 
understanding  of  data  analysis  would  provide  an 
opportunity  for  improved  accuracy,  error  estimation,  and 
data  reduction  rates. 


THEORY  OF  THRUST  STAND  DYNAMICS 

The  analytical  model  of  the  thrust  stand  is 
chosen  to  be  that  of  an  underdamped  oscillator.  The 
general  equation  governing  the  oscillatory  motion  of  the 
thrust  stand  under  load  is  that  of  a  second  order  damped 
harmonic  oscillator  with  constant  coefficients,  given  in 
Eq.  1.  The  displacement  (x)  is  a  function  of  time  (t), 
damping  rate  (a=2->/ Kf  m  oc ),  spring  constant  (K), 
effective  system  mass  (m),  thruster  load  ( F ),  and  initial 
conditions  for  displacement  and  velocity. 

m  m 

The  modeling  effort  first  examines  thrust  stand 
response  under  uniform,  periodic,  impulsive  load.  This 
case  is  representative  of  an  ideal  system  where  neither  the 
thruster,  nor  the  thrust  stand  exhibits  any  non-uniform 
characteristics;  the  displacement  data  is  accurate  and  the 
thruster  provides  a  reproducible  impulse  bit. 
Representative  displacement  oscillations  are  shown  in 
Figure  1  for  three  values  of  x.  Equations  describing  thrust 
stand  displacement  as  a  function  of  time  and  operating 
parameters  are  presented  below.  Next,  the  effect  of  noisy 
displacement  data  on  impulse  bit  calculations  is 
examined.  Finally,  the  relationship  between  measured 
displacement  and  thrust  for  non-uniform,  periodic, 
impulsive  loads  is  considered. 

With  the  assumption  of  a  periodic  impulsive 

load,  the  oscillatory  motion  is  dependent  upon  a,  x,  K,  the 
period  between  impulses  (tp),  and  the  thruster  impulse  bit 
(fbit).  The  equations  governing  thrust  stand  response  then 
reduce  to  those  given  in  Eqs.  2-4.  The  fractional  time 
within  an  impulse,/  defines  the  time-wise  location  within 
a  given  impulse.  The  pulse  period  fraction,  / ,  is  equal  to 
fraction  of  time  within  the  impulse  that  occurs  after  the 
displacement  maximum  (Fig.  1).  /  is  calculated  using 
Eqs.  2  and  3;  the  displacement  at  /  is  equal  to  xmin.  Eq.  4 
is  used  to  compute  the  maximum  thrust  stand 
displacement  Cw).  Eq.  2  is  used  to  calculate  the  time 
history  of  displacement,  which  begins  with  an  impulse  at 

/  =  /-!,  and  ends  at  /  =  / . 


r(/.T,a)  =  *max* 


-2afKT 


(2) 


.v(/,r,a)=x(/-l,T,a) 


(3) 
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Fig.  1  Idealized  thust  stand  displacement 
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Sin^27t f^J\-a2 rj<?  2nfar  + 
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(4) 


From  the  results  obtained  above,  the 
displacement  associated  with  the  average  thrust  can  be 
computed  (Eq.  5).  Eqs.  4  and  5  show  that  the  average 

displacement,  xmr,  is  equal  to  this  wil1  later  be 

shown  to  be  a  special  case  of  oscillator  behavior. 


A.  Displacement  Noise 

The  effect  of  displacement  data  noise  on  impulse 
bit  calculations  is  examined  through  the  addition  of 
random  noise  to  the  ideal  displacement  curves  described 
by  Eqs.  2-4.  Displacement  noise  is  computed  using  a 
random  amplitude  with  standard  deviation  equal  to  a 


fractional  error,  y,  multiplied  by  the  amplitude  of 
maximum  displacement.  The  resulting  displacement  data 
is  analyzed  using  a  commercial  non-linear  Fitting 
algorithm10.  A  modified  version  of  Eq.  2,  which  allows 
for  shifts  in  time  of  maximum  displacement,  is  the  fitted 
to  the  data.  The  algorithm  Fits  four  independent 
parameters  and  computes  the  oc,  x,  f ,  xmax,  then  for 
each  impulse. 

B.  Non-Uniform  Impulse  Loads 

The  effect  of  non-uniform,  periodic,  impulsive 
loads  on  average  thrust  calculations  is  examined.  This 
analysis  assumes  that  disturbance  forces,  the  effect  of 
noise,  and  drift  associated  with  thrust  stand  measurement 
are  negligible.  A  transient  solution  to  Eq.  1  was 
implemented,  with  the  impulse  bit  for  each  pulse 
computed  using  a  random  addition  to  the  average  with  a 
standard  deviation  equal  to  a  fractional  variation,  8,  from 
the  average  (Fig.  2).  An  analytic  solution  for  displacement 
under  non-uniform,  periodic,  impulsive  load  was  derived. 
The  general  solution  for  damped  oscillator  motion  is 
given  in  Eqs.  6-8.  The  recursion  relation  for  the  fitting 
parameters  (Eqs.  9-10)  was  derived  for  an  underdamped 
thrust  stand.  The  parameter  /  again  describes  the 
dimensionless  time  within  an  impulse  cycle  of  duration  tp. 
For  this  case, /varies  from  0  to  1.  Impulses  are  numbered 
by  the  index  J  in  parenthesis.  The  first  derivative  term  in 
Eq.  10  is  with  respect  to  the  dimensionless  parameter/. 

.v(q  ,c2 ,/)  =  c,  yj  (a,  T,/)  +  c2  y2  (a,T,/)  (6) 


n  {a,x,f)=e~2arKf  Cos^  l-a2/^rj  (7) 

y2  {a,T,f)=e~2atKf  Sin^\-a2  fxr  j  (8) 
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Fig.  2  Simulated  thrust  stand  displacement  with 
impulse  bit  variance  (5  =  0.05) 

Eqs.  6-10  were  examined  analytically  to 
determine  the  relationship  between  the  average 
displacement  and  the  average  impulse  bit 

t  iave=- L  y  //.  )  for  an  impulse  set  of  N  impulses.  The 
(  bit  N  flxblt 

recursion  relation  is  solved  using  a  summation  and  matrix 
inversion  technique.  It  is  found  that  the  impulse  induced 
displacement  uncertainty  (£^),  the  fractional  difference 

between  the  measured  average  displacement  and  the 
displacement  associated  with  the  average  thrust 
(/'^  / 1  ),  is  equal  to  Eq.  12;  equation  simplification 

was  accomplished  using  commercial  symbolic 
mathematical  manipulation  routines1".  The  parameter  dt 
and  d2  are  defined  in  Eqs.  13-15.  dk  equals  the  difference 
between  the  beginning  and  ending  coefficients  divided  by 
C^eal,  the  value  for  uniform  impulses  (8=0).  Thus,  when 
the  starting  and  ending  displacement  curves  are  equal,  as 
is  the  case  for  uniform  periodic  impulses,  the  average 
displacement  is  equal  to  the  displacement  associated  with 
the  average  thrust. 


,ave 

‘m. 

Ktn 
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,ave 
' 'bit 
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This  analysis  proves  that,  subject  to  the  model 
assumptions,  average  displacement  measurements 
approximate  the  displacement  associated  with  average 
thrust  to  an  accuracy  that  is  inversely  proportional  to  the 
number  of  impulse  bits  averaged.  Also,  £d  is  equal  to 
zero  for  any  set  of  impulses  for  which  the  first  and  last 
displacement  histories  are  the  same. 
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NUMERICAL  ANALYSIS 


The  displacement  ( 


*  bit  s 
Ktp 


associated  with  the 


average  thrust  relative  to  xmax  and  xmj„  is  represented  by 
the  thrust  coefficient  ( (j) )  (Eq.  16).  The  values  for  xnmx 


and  xmi„  are  calculated  using  Eq.  2  for  /  =  0  and  /  -  /, 
respectively.  By  characterizing  (p  for  a  range  of  damping 
coefficients  and  pulse  period  ratios,  the  displacement 
associated  with  average  thrust  can  be  computed 
accurately  with  the  measurement  of  x„u„,  xmi„,  a,  and  x. 
Using  Eqs.  2-5,  <j)  was  calculated  as  a  function  of  a  and  x 
(Fig.  3).  The  figure  shows  that  the  thrust  coefficient  varies 
only  slightly  for  pulse  period  ratios  less  than  0.5. 


<p(a.x) 


~-vmin 
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Fig.  3  Thrust  coefficient  as  a  function  of  a 
and  X 


For  the  special  case  of  undamped  oscillation,  a 
simple  expression  exists  for  <p  (Eq.  17).  When  this 


analytic  expression  is  used 


to  compute 


I  bit 
Ktp ’ 


the 


fractional  error  in  impulse  bit  resulting  from  the 
approximation  is  less  than  2%  for  pulse  periods 
approaching  0.4  (Fig.  4). 


<p{  0,t) 
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Fig.  4  Fractional  thrust  error  with 
assumption  of  no  damping 


An  improved  approximation  for  <j)  was  determined 
using  Eq.  17  as  the  first  term  in  a  functional  expansion. 
Calculated  values  of  <p  were  used  to  determine  a 
polynomial  fit  to  the  expression.  A  36  by  36  matrix  of 
calculated  values  for  a  (0.1  -0.9)  and  x  (0.01-0.7)  provided 
the  input  to  the  regression  algorithm10.  The  fitted  function 
includes  first  and  second  order  terms.  The  resulting 
expression  is  given  below: 

<p(a.r)  =  0(0, t)  + 

(-0.0098881  +  a(0. 014532)  +  r(0.1 1039)  +  (18) 

a  r  (-0.1 7889)  +  a2  (0.0042004)  +  x2  (-0.14525) 

This  expression  provides  a  significant  improvement  in 
accuracy  over  the  strictly  analytic  approximation.  Use  of 
this  approximation  provides  a  thrust  error  less  than  2%  for 
all  values  of  a  and  x  shown  in  Fig.  5. 
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Fig.  5  Fractional  thrust  error  with 
second  order  polynomial  approximation 


around  5%  for  x  =  0.3.  These  results  are  consistent  with 
the  increase  in  the  effect  of  damping  at  large  values  of 
x,  and  the  anticipated  sensitivity  to  noise  associated  with 
small  values  of  X  where  damping  is  less  of  a  perturbing 
force.  No  significant  influence  of  X  on  the  other  fitting 
variables  was  observed.  For  all  cases  examined,  the 
average  value  for  the  differences  was  at  least  two  orders 
of  magnitude  less  than  the  standard  deviation,  indicating 
that  the  effect  of  random  noise  is  distributed  evenly  about 
the  exact  value. 


Ey  =  Std  dev | 


J\calc  _j  exact 


bit 


bit 
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bit 


(19) 


The  results  show  that  er  decreases  as  the  pulse 
period  increases  within  the  range  from  0.1  to  0.3.  This 
effect  is  consistent  with  the  large  variability  in  calculated 
values  of  a  for  small  values  of  X.  Within  the  region 
between  x  =  0.3  and  0.7,  where  a  is  less  variable,  £y  is 
approximately  constant. 


ERROR  ANALYSIS  RESULTS 

This  paper  examines  the  thrust  measurement 
uncertainties  associated  with  two  types  of  error  sources, 
random  noise  and  non-uniform  impulse  bit  sets.  To 
evaluate  the  magnitudes  associated  with  these  error 
sources,  two  different  analysis  methods  were  used. 
Displacement  uncertainty  associated  with  noisy  thrust 
stand  displacement  data  was  analyzed  by  fitting  simulated 
random  thrust  stand  displacement  perturbations.  The 
effect  of  non-uniform  impulse  bit  thruster  operation  is 
evaluated  using  numerical  simulation  and  the  analytic 
expression  for  thrust  measurement  error  given  in  Eq.  12. 


A.  Displacement  Noise 

Displacement  noise  perturbation  magnitude  was 

varied  from  zero  to  4.0x10  times  the  peak  oscillation 
amplitude  for  each  case.  The  factor  y  is  the  fractional 
magnitude  (standard  deviation)  of  the  noise  with  respect 
to  peak  oscillation  amplitude.  Results  from  the  fitting  of 
50  impulses  for  each  condition  are  shown  in  Fig.  6.  For 
all  cases,  a  damping  coefficient  of  0.3  was  used.  The 
displacement  uncertainty  for  noisy  displacement  data,  £r 
is  defined  as  the  standard  deviation  of  the  set  of  fractional 
differences  between  the  fitted  and  exact  impulse  bits  (Eq. 
19).  The  fitted  impulse  bit  is  computed  using  Eq.  5  and 
the  fitted  values  for  a,  x,  /,  and  xmax.  The  fitting  of  the 
damping  coefficient  was  seen  to  depend  strongly  on 
X,  with  standard  deviation  variations  of  more  than  100% 
above  the  average  resulting  for  x  =  0.1  compared  to 


Fig.  6  Displacement  uncertainty  for  noisy 
displacement  data  (oc  =  0.3) 

In  addition  to  the  a  =  0.3  results,  additional  calculations 
were  made  for  a  =  0.7,  x  =  0.1.  At  this  high  level  of 
damping,  the  larger  variability  of  calculated  values  of  oc 
was  also  observed  for  x  =  0.1.  Though  the  calculated 
value  of  £y  did  not  differ  significantly  from  the  results 
shown  in  Fig.  6,  the  number  of  iterations  required  to 
produce  a  converged  solution  increased  significantly  from 
that  observed  for  ct  =  0.3  (from  about  10  to  over  100). 
This  is  consistent  with  the  large  variability  in  the  fitted 
value  of  oc,  which  results  in  values  of  ct  that  exceed  unity, 
the  range  of  applicability  of  the  equations  used  to  model 
thrust  stand  response. 
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B.  Non-Uniform  Impulse  Loads 

The  effect  of  non-uniform  impulse  bits  on  thrust 
measurement  was  evaluated  numerically  using  Eq.  12  and 
numerical  simulations  of  thrust  stand  response.  The 
measurement  uncertainty,  £,, ,  associated  with  randomly 
variable  impulse  bit  is  calculated  from  the  standard 
deviation  of  a  set  of  Ed  values.  The  parameter  d,  and  d2 
are  calculated  for  each  impulse;  a  first  order  examination 
of  the  distribution  of  these  parameters  confirmed  that  they 
are  equally  distributed  about  dk  =  0.  For  the  results 
presented  in  Figs.  7  and  8,  each  point  was  the  result  of 
100  data  sets.  The  dependency  of  Ed  on  a  and  X  is 
presented  in  Fig.  7.  For  these  cases,  the  impulse  sets 
consisted  of  5  impulses  (N=5)  and  the  impulse  bit 
standard  deviation  (5)  is  5%.  These  results  show  that,  in 
general,  Ej  decreases  with  increasing  a  and  x;  though  the 
a  =  0.7  data  conflicts  somewhat  with  this  conclusion,  the 
limited  sampling  sizes  used  can  explain  such  a 
discrepancy  (discussed  below).  The  largest  variation 
occurs  for  a  and  x  between  0. 1  and  0.3. 


impulses  per  pulse  set,  100  sets) 

The  effect  of  impulse  set  length  and  impulse  bit 
variability  is  shown  in  Fig.  8.  As  expected  from  Eq.  12, 
the  results  show  a  near  inverse  relation  between  £,/  and  N. 
In  addition,  £,/  is  approximately  proportional  to  impulse 
variability  ( S ).  These  results  show  that  thrust 
measurements  using  integrated  displacement  data  can 
accurately  represent  the  average  thrust  for  even  relatively 
small  impulse  set  sizes.  For  example,  the  integration  of  10 


impulses  is  sufficient  to  provide  a  less  than  2% 
measurement  standard  deviation  for  a  random  10% 
impulse  bit  standard  deviation. 

Based  on  calculations  with  large  sampler  sizes 
(200,  400,  and  800)  compared  to  the  100  used  for  the 
results  shown  in  Fig.  7  and  8,  the  accuracy  of  the  data 
presented  is  estimated  to  be  within  50%.  An  accuracy  of 
that  magnitude  is  consistent  with  the  relative  location  of 
the  tx  =  0.7  line  in  Fig.  7. 


Fig.  8  Displacement  measurement 
uncertainty  for  non-uniform  periodic  impulse 

bits  (a  =  0.3,x  =  0.1) 

DISCUSSION 

The  analytic  examination  of  thrust  stand 
dynamics  provides  sensitivity  estimates  for  thrust  stand 
measurement  accuracy  as  a  function  of  data  reduction 
method.  The  results  presented  show  the  relative 
sensitivity  of  thrust  stand  accuracy  on  damping  and  the 
ratio  between  pulse  period  and  thrust  stand  natural 
(undamped)  period. 

The  data  show  that  the  uncertainty  associated 
with  random  displacement  noise  decreases  as  the  pulse 
period  ratio,  x,  is  increased  when  a  curve  fitting  technique 
is  used  to  calculate  individual  impulse  bits,  this  effect  is 
greatest  for  x<0.3.  This  technique  also  experiences 
convergence  difficulties  when  a  approaches  unity. 

An  analytic  expression  for  measurement 
accuracy  was  derived  for  an  underdamped  thrust  stand 
under  non-uniform  period  impulsive  loading.  This 
expression  provides  a  means  of  estimating  the  uncertainty 
associated  with  a  very  simple  method  of  thrust 
measurement  analysis,  the  averaging  of  thrust  stand 
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displacement.  This  analysis  indicates  that  measurement 
uncertainties  associated  with  relatively  large  impulse  bit 
variability  can  be  reduced  through  the  integration  of  a 
relatively  small  number  of  impulses,  on  the  order  of  10. 
With  such  a  small  number  of  impulses  required  to 
accurately  determine  average  thrust,  it  is  believed  that  this 
method  would  be  applicable  for  the  analysis  of  even 
transient  effects,  such  as  the  rapid  increase  in  thruster 
impulse  bits  due  to  discharge  energy  change.  In  addition, 
data  analysis  is  expected  to  be  very  rapid  due  to  simplicity 
of  the  numerical  technique;  the  displacement  data  is 
averaged  with  respect  to  time.  The  data  also  show  that 
measurement  uncertainty  decreases  as  X  is  increased;  this 

effect  is  greatest  for  x  <  0.3. 

The  approximations  used  in  the  displacement 
averaging  thrust  measurement  technique  are  believed  to 
accurately  represent  the  thrust  stand  system.  While 
random  and  periodic  disturbances  are  omitted  from  the 
analysis,  it  is  anticipated  that  they  will  largely  cancel 
provided  that  their  period  of  oscillation  is  much  less  than 
the  averaging  period.  Displacement  drift  is  not  expected 
to  be  a  significant  source  for  measurement  uncertainty 
provided  it  remain  linear  over  the  duration  required  for 
displacement  averaging. 

Two  issues  remain  that  will  require  experimental 
confirmation  to  determine  their  significance.  The  first  is 
non-random  electromagnetic  induced  displacement  noise. 
Effective  shielding  of  data  acquisition  electronics  is 
expected  to  adequately  address  this  problem.  The  second 
issue  is  that  of  the  assumption  of  constant  coefficients  in 
the  second  order  differential  equation.  While  not  expected 
to  be  an  issue,  the  derivation  is  only  exact  if  the 
coefficients  are  constant.  It  is  anticipated  that  the 
accuracy  of  the  technique  can  be  confirmed  through  the 
comparison  of  thrust  analysis  results  for  a  range  of  pulse 
period  ratios.  Variation  of  both  thruster  pulse  period  and 
thrust  stand  natural  period  would  be  used. 

An  interesting  conclusion  that  can  be  drawn  from 
this  analysis  is  that  there  is  an  accuracy  advantage 
associated  with  the  use  short  period  (x  >  0.3)  thrust  stand 
designs.  With  T  proportional  to  Kim  ,  this  indicates  that 
more  accurate  measurements  can  be  obtained  using  less 
massive  (also  more  compact,  since  m  is  proportional  to 
the  thrust  stand  moment  of  inertia6)  thrust  stands.  Such 
systems  could  be  used  in  smaller  vacuum  chambers  and, 
due  to  their  reduced  size,  would  be  simpler  to 
mechanically  isolate.  A  small,  well  isolated,  high 
accuracy  thrust  stand  could  provide  a  means  of  accurately 
measuring  impulse  bits  in  the  sub-lOOpN-s  range. 


SUMMARY 

This  paper  presents  a  numerical  analysis  of 
thrust  stand  response  during  pulsed  thruster  operation.  A 
summary  of  the  results  is  given  below; 


•  Displacement  oscillation  shape  as  a  function  of  a  (0. 1 
to  0.9),  and  x  (0.05  to  1.5)  are  presented.  Impulse 
approximation  techniques  accurate  to  within  2%  for  x  less 
than  0.7  are  shown. 

•  Displacement  signal  noise  effects  were  quantified  for 
(x=0.3  with  x  varying  from  0.1  to  0.7.  Thrust  measurement 
uncertainties  below  2%  result  from  fractional  noise 
amplitudes  less  than  0.004.  Measurement  uncertainty  due 
to  noise  decreases  when  x  is  increased  to  -  0.3. 

•  An  analytic  expression  relating  average  thrust  and 
integrated  thrust  stand  displacement  was  derived. 
Numerical  simulations  were  used  to  estimate  thrust 
uncertainty  dependence  on  the  number  of  impulses  used 
for  averaging.  Measurement  uncertainty  due  to  impulse 
bit  variability  decreases  when  x  is  increased  to  ~  0.3. 
Integration  of  10  impulses  is  adequate  to  achieve 
measurement  uncertainty  below  2%  for  a  random  10% 
impulse  bit  variation. 
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